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Pion-pion and pion-proton correlations - new results from CERES* 

Dariusz Antonczyk and Dariusz Miskowiec for the CERES Collaboration 
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Results of a new two-particle correlation analysis of central Pb+Au collision data at 158 GeV per nucleon are 
presented. The emphasis is put on pion-proton correlations and on the dependence of the two-pion correlation 
radii on the azimuthal emission angle with respect to the reaction plane. 
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I. EXPERIMENT AND DATA ANALYSIS 

CERES is a dilepton experiment at the CERN SPS, known 
for its observation of enhanced production of low mass e + e~ 
pairs in collisions between heavy nuclei [ 1 ] . The upgrade of 
CERES in 1997-1998 by a radial Time Projection Chamber 
(TPC) allowed to improve the momentum resolution and the 
particle identification capability while retaining the cylindri- 
cal symmetry (Fig. [TJ. The TPC also opened the possibility 
of measuring hadrons. The upgraded experiment has been de- 
scribed in detail elsewhere y2|]. 

The measurement of central Pb+Au collisions at the maxi- 
mum SPS energy of 158 GeV per nucleon in the fall of 2000 
was the first run of the fully upgraded CERES and at the same 
time the last run of this experiment. About 30 million Pb+Au 
collision events at 158 GeV per nucleon were collected, most 
of them with centrality within the top 7% of the geometrical 
cross section Og = 6.94 b. Small samples of the 20% and the 
minimum bias collisions, as well as a short run at 80 AGeV, 
were recorded in addition. The first two-particle correlation 
analysis performed on a subset of these data resulted, among 
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FIG. 1 : Upgraded CERES setup in 2000. The apparatus has a cylin- 
drical symmetry. The beam enters from the left. The silicon detectors 
(SDD) give tracking and vertex reconstruction, and the Ring Imag- 
ing Cherenkov detectors (RICH) electron identification. The results 
presented here are mostly based on the momentum (and energy loss) 
measurements performed with the Time Projection Chamber (TPC). 



others, in an improved procedure to account for the Coulomb 
interaction yfl and a new postulate of an universal freeze-out 
criterion [4]. The current analysis features a better momentum 
resolution 



Ap 



= 2%® 1%-p/ (GeV/c), 



(1) 



a better understanding of the two-track resolution, and was 
performed on the full data set. The n(p t ,y) distribution of 
the analyzed pairs is shown in Fig. [2] The two-track reso- 
lution cuts applied to the true pairs and to the pairs from event 
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FTG. 2: Distributions of the analyzed two-pion (top) and pion-proton 
(bottom) pairs. The midrapidity is ygl2=2.9\. 
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FIG. 3: Track reconstruction efficiency in the case of the magnetic 
field bringing the tracks apart from (left) or closer to (right) each 
other. The two topologies were dubbed "sailor" and "cowboy", re- 
spectively, and required different two-track separation cuts A0: 38- 
45 mr for sailor and 90-140 mr for cowboy, depending on the trans- 
verse momentum. The separation cut was 8-9 mr. 



mixing were different for the two possible track pair topolo- 
gies (Fig. [3). It should be noted that the required two-track 
cuts depended somewhat on the quality cuts applied to single 
tracks: the higher number of hits required for single tracks, 
the more pairs were lost because of the finite two-track res- 
olution. The two-pion analysis was performed in the longi- 
tudinally co-moving frame (LCMS) defined by the vanishing 
z-component of the pair momentum. The momentum differ- 
ence in this frame, q = P2-P1, was decomposed into the "out", 
"side", and "long" components following the Bertsch-Pratt 
convention, with q out pointing along the pair transverse mo- 
mentum and <7i ong along the beam. The non-identical correla- 
tions were analyzed in the pair c.m.s., the frame in which the 
pair momentum is zero. The two components there, and 
q±, were defined such that was along the pair transverse 
momentum, i.e. was equal to gout if the latter is calculated 
in the pair c.m.s. 



II. TWO-PION CORRELATIONS 



F c (<7inv) factor, qx m = J — — p^) 2 , accounts for the mutual 
Coulomb interaction between the pions and was calculated by 
averaging the nonrelativistic Coulomb wave function squared 
over a realistic source size. The Coulomb factor was attenu- 
ated by X similarly as the rest of the correlation function peak; 
the importance of this approach was demonstrated in J^. The 
fits were performed by the minimum negative loglikelihood 
method with the Poissonian number of true pairs and were 
done separately for each pair (p t ,y) bin. The HBT radii ob- 
tained from the fit were corrected for the finite momentum 
resolution. The correction was determined by Monte Carlo 
and was rather insignificant for R s ^ e and Ri ong ; for R out it gets 
as large as sa 20% for the highest bin of the pair p t . The ob- 
tained HBT radii show a strong p t dependence (Fig. [4]). The 
^side and ^iong radii were fitted with [6, 7] 
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with the freeze-out temperature T fixed to be 120 MeV. The 
results of the fit are shown in Fig. [5] The obtained expansion 
time x of 5.5-6.5 fm, geometrical source size Rq of 6-8 fm, 
and transverse expansion rapidity r\f of 0.7-0.8, in reasonable 
agreement with the results of the previous analysis of a subset 
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The 7T.~7r,~ and n + n + correlation functions, defined as the 
three-dimensional distributions of pion pairs from the same 
event n(q), normalized to the analogous distributions of pairs 
constructed from different events (event mixing), were fitted 
by 



C 2 (q)=N-i{l-X) + X-F c (q im ) 



1+exp^- £ Rijqiq/J 
(2) 

The normalization factor N is needed because the number of 
pairs from event mixing is arbitrary. The correlation strength 
X < 1 reflects the tails of the source distribution caused by the 
pions from long-lived resonances, the finite q-resolution, and 
the contamination of the pion sample by other particle species. 
The Rfj fit parameters, with the indices ij being {out, side, 
long}, are related to the size of the source emitting pions of 
given momentum [5] and will be called here HBT radii. The 
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FIG. 4: Transverse momentum dependence of the source radii ob- 
tained from the ir~Tt~ and Jl + Jt + correlation analysis. The top and 
the bottom panels show the most and the least central bins (0-2.5% 
and 15-35% of Oq), respectively. The shown in the abscissa is 
the pion pair transverse momentum divided by two. The data are 
preliminary and come from central Pb+Au collisions at 158 GeV per 
nucleon. 
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FIG. 5: Centrality dependence of the pion source parameters obtained from fits ((3) and ©. The data are preliminary and come from central 
Pb+Au collisions at 158 GeV per nucleon. 



of data 1 3], indicate a long-living, longitudinally and transver- 
sally expanding pion source. 

III. AZIMUTHAL ANGLE DEPENDENT HBT 

A fireball created in a collision with a finite impact param- 
eter is elongated in the direction perpendicular to the reaction 
plane. In the course of expansion, with the pressure gradient 
larger in-plane than out-of-plane, the initial asymmetry should 
get reduced or even reversed. A dependence of the two-pion 
correlations on the pair emission angle with respect to the re- 
action plane would be a signature of the source eccentricity 
at the decoupling time. The azimuthal angle of the reaction 
plane was estimated by the preferred direction of the particle 
emission aka elliptic flow. The particles were weighted with 
their transverse momentum: 



Qi = I>cos(2^) 

i 

Qi = Y*P' sin ( 2 0i)- 



(5) 
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The raw distribution «(2f , Q2 ) i s shown in Fig. [6] The reac- 
tion plane angle was calculated (modulo Jt) from the calibrated 
Q2 components via 



<I>rp = - arctan 



Ql 
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The resolution of the so determined reaction plane angle, es- 
timated via the subevevent method, was 31-38°. 

With the event plane known (within the resolution) event- 
by-event the pion pairs can be sorted into 8 bins covering 
(— 7t/2,7t/2) according to their azimuthal angle with respect 
to the reaction plane <t>* = <I>p a i r — ^rp. During event mix- 
ing it was required that the two events had a similar reaction 
plane angle ^Rp. The eight correlation functions were fitted 
as described in Section lU and the resulting R out , 7? s ide, ^long, 



and the cross-terms are plotted versus <f>*. The squared source 
radii were then fitted with 



Rj=Rl + 2Rl 2 cos(2<P* 



(8) 



with i denoting {out,side,long}. While the /?,,o's obtained co- 
incide with the results of the standard HBT analysis presented 
in Section HI1 the second Fourier components Ri z's represent 
the eccentricity of the observed pion source. The normalized 
second Fourier component of R s id e is shown in Fig. [7] As far 
as the limited centrality range allows to judge the measured 
R s id & anisotropy is consistent with zero, in contrast to the out- 
of-plane elongated pion source observed both at the AGS lH|] 
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FIG. 6: The raw distribution of the (0^,0^) vector used to deter- 
mine the azimuthal orientation of the reaction plane. The distribution 
was recentered and made round run-by-run in order to make the re- 
sulting distribution of the reaction plane angle uniform. 
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and at RHIC & The source eccentricity, thus, seems to be 
joining the exclusive club of heavy ion observables which be- 
have non-mono tonic ally with the collision energy, the other 
members being the flow and the strangeness-to-entropy ratio. 
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FIG. 7: Azimuthal source eccentricity, represented by the normal- 
ized second Fourier component of ^ s ide('I ) *)> vs. centrality. The 
CERES result (black full dots) is much closer to zero than the anal- 
ogous measurement at the AGS (open red squares) and RHIC (open 
blue triangles). 



IV. PION-PROTON CORRELATIONS 

For the top central 7% of the geometrical cross section the 
high event statistics allows to perform the proton-pion corre- 
lation analysis. The shape of non-identical particle correlation 
functions C(q) reflects the shape of the relative source distri- 
butions S(r^ — /j'). Particularly, a difference between the aver- 
age freeze-out position or time of two particle species reveals 
itself as an asymmetry of the correlation function at small q 
ifioll . A two-dimensional 7T,~p correlation function C{q\\,q±) 
and its slice C(q\\) for q±<50 MeV/c are shown in Fig. [8] 
The peak at low q comes from the attractive Coulomb inter- 
action. The peak asymmetry indicates that the proton source 
is located at a larger radius than the pion source, or that pro- 
tons are emitted earlier than pions. The asymmetry can be 
conveniently parametrized by fitting a Lorentz curve, which 
happens to match the shape, separately to the left and to the 
right half of the peak, and taking the ratio of the two widths: 
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where AHs a normalization factor and a is the peak amplitude. 
The asymmetry is then defined as the ratio between the two 
widths A = G_/o"+. 
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FIG. 8: Two-dimensional 7t~p correlation function C(qpq±) (left) 
and its slice C(q^) for q^<50 MeV/c (right). The peak asymmetry 
in the right hand plot indicates that the proton source is located at a 
larger radius than the pion source, or that protons are emitted earlier 
than pions. The asymmetry can be parametrized using Eq. (9}. 



The asymmetry was translated to a spatial displacement be- 
tween the proton and pion sources using a Monte Carlo pah- 
generator with realistic source sizes and the Coulomb wave 
function squared as a weighting factor for each pair. The re- 
lation between the two quantities is shown in Fig. [9] The re- 
sulting displacement Ax as a function of the transverse pair 
momentum is shown in the right panel of Fig. |T0j In the left 
panel the raw asymmetry A is shown. The asymmetry and the 
displacement derived from it vanish in the limit of small pair 
p t as is expected for symmetry reasons. The displacement is 
rather similar to the spatial displacement seen in UrQMD v. 1 .3 
111] (green line). The red and blue lines represent a simulta- 
neous fit to R s itie(Pt) and Ax(p t ) using, respectively, Eq.fO and 
the formula derived in 11211 : 
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FIG. 9: Relation between the pion-proton correlation peak asymme- 
try and the displacement between the source of protons and pions, 
obtained from Monte Carlo. This relation was used to translate the 
measured asymmetries to displacements. 
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FIG. 10: Left: 7t + -proton (top) and Jt~-proton (center) correlation peak asymmetry, defined as the ratio between the two half-widths in Eq. 
(9), vs. the pair transverse momentum. The asymmetry is larger than one which indicates that protons on average freeze-out at a larger radius 
(or earlier time) than pions. No asymmetry is visible in the Jt + Jt~ correlations (bottom panel). For symmetry reasons all peak asymmetries 
must vanish at small p,'s. Right: the source displacement deduced from the asymmetry. The green lines represents the displacements extracted 
from UrQMD for p t 's between and 2 GeV/c. The long lines show the fit described in text. The data are preliminary and come from central 
Pb+Au collisions at 158 GeV per nucleon. 



with the pair transverse mass 



m t = 



i 



m p + m K \ 2 
and the pair transverse velocity 



V p i- J 



(12) 



m K + m p \ 2 

(11) The relation between the rapidity r|^ and the velocity po of 
transverse flow, which occur in the R s ^ s and Ax fit formulas, 
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FIG. 11: Transverse momentum spectra of negatively charged particles, K , and 0-meson (top) and elliptic flow coefficient for pions, K , and 
A (bottom), compared to the blast-wave model (red line). The two green lines in the top left plot represent negative pions and kaons, the sum 
of the two being the red line to be compared to the data. The data are preliminary and come from central Pb+Au collisions at 158 GeV per 
nucleon.The black and blue points represent the negative and positive pions, respectively. 



respectively, is 



1, 1+Po 



(13) 



The fit yields a p of 0.65-0.70 and a R G of about 7.5 fm. The 
fit is dominated by the R s a e data with their small error bars 
but it reproduces the Ax values quite well. This indicates that 
the finite displacement between the sources of pions and pro- 
tons has a similar origin as the p, dependence of /? s ide> namely 
the transverse flow. Other, possibly more interesting, contri- 
butions to it cannot be addressed with the present statistics of 
the data. 



good. Two of the eight parameters of the model were fixed: 
the freeze-out temperature T = 100 MeV and the sharpness of 
the emission volume a = 0.01 (sharp). The other six param- 
eters, adjusted to the data by the simplex method, include the 
transverse flow rapidity p of 0.88 and 0.85 in-plane and out- 
of-plane, respectively; the source radius (R X ,R } ) of (11.26, 
11.42) fm; the emission time (i.e. the duration of the expan- 
sion) x of 7.4 fm, and the duration of emission Ax = 1 .55 fm/c. 
The nearly-circular source R x « R y was enforced by the small- 
ness of the second Fourier component of the HBT radii (bot- 
tom part of Fig. [T2l . 



VI. SUMMARY 



V. BLAST- WAVE PARAMETRIZATION 

While it is notoriously difficult to describe many ob- 
servables of heavy ion collisions within the same theoreti- 
cal approach, in recent years simple hydrodynamics-inpired 
parametrizations turn out to be quite successful in this as- 
pect. The blast-wave model of 01311 . with its 8 parameters 
adjusted accordingly, nicely reproduces the transverse mo- 
mentum spectra, the elliptic flow, and the two-particle corre- 
lations of CERES, including the emission angle dependence 
(or, rather, the lack of it). The transverse spectra and the ellip- 
tic flow coefficients from the top 7% of Og for several hadron 
species are shown along with the blast-wave model lines in 
Fig-El in Fig.[T2lthe two-particle correlation results are com- 
pared to the calculations performed with the blast-wave model 
with the same parameter values. The agreement is rather 



In summary, the new high-statistics CERES hadron data 
allow for state of the art analysis of two-particle correlation 
data. The pion source anisotropy, accessible via the depen- 
dence of the two-pion HBT radii on the emission angle with 
respect to the reaction plane, is unexpectedly small compared 
to the analogous results obtained at lower and higher ener- 
gies. The pion-proton correlations indicate that the proton 
source is located at a larger radius than the pion source, or 
that protons are emitted earlier than pions. The amount of 
the displacement, and its functional dependence on the trans- 
verse momentum of the two involved particles, corroborate 
the transverse expansion picture deduced from the behavior 
of the R s ide HBT radius. Finally, the blast-wave parametriza- 
tion is able to describe simultaneously the transverse spectra 
and the elliptic flow coefficients of several hadron species, and 
the two-particle correlations. 
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FIG. 12: CERES HBT radii squared and the pion-proton displacement (top), HBT cross-terms (middle), and the second Fourier components 
of the HBT radii vs. emission angle with respect to the reaction plane (bottom), compared to the blast-wave model calculations (red line). The 
data are preliminary and come from central Pb+Au collisions at 158 GeV per nucleon.The black and blue points represent the negative and 
positive pions, respectively. 
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